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Abstract. We introduce a new invariant, the coronal of a graph, and use 
it to compute the spectrum of the corona G o H of two graphs G and H. 
In particular, we show that this spectrum is completely determined by the 
spectra of G and H and the coronal of H. Previous work has computed the 
spectrum of a corona only in the case that H is regular. We then explicitly 
compute the coronals for several families of graphs, including regular graphs, 
complete n-partite graphs, and paths. Finally, we use the corona construction 
to generate many infinite families of pairs of cospectral graphs. 



1. Introduction 

Let G and H be (finite, simple, non-empty) graphs. The corona G o H oi G 
and H is constructed as follows: Choose a labeling of the vertices of G with labels 
1,2, ... ,m. Take one copy of G and m disjoint copies of H, labeled Hi, ... , H„i, 
and connect each vertex of Hi to vertex i of G. This construction was introduced 
by Frucht and Harary [5] with the (achieved) goal of constructing a graph whose 
automorphism group is the wreath product of the two component automorphism 
groups. Since then, a variety of papers have appeared investigating a wide range 
of graph-theoretic properties of coronas, such as the bandwidth [2], the minimum 
sum [12], applications to Ramsey theory |3, etc. Further, the spectral properties 
of coronas are significant in the study of invertible graphs. Briefly, a graph G is 
invertible if the inverse of the graph's adjacency matrix is diagonally similar to 
the adjacency matrix of another graph G"*". Motivated by applications to quan- 
tum chemistry, Codsil [4] studies invertible bipartite graphs with a unique perfect 
matching. In response to his question asking for a characterization of such graphs 
with the additional property that G = G"*", Simian and Cao |10| determine the 
answer to be exactly the coronas of bipartite graphs with the single-vertex graph 
Ki. 

The study of spectral properties of coronas was continued by Barik, et. al. in [T], 
who found the spectrum of the corona G o if in the special case that H is regular. 
In Section 2, we drop the regularity condition on H and compute the spectrum of 
the corona of any pair of graphs using a new graph invariant called the coronal. In 
Section 3, we compute the coronal for several families of graphs, including regular 
graphs examined in [T], complete n-partite graphs, and path graphs. Finally, in 
Section 4, we see that properties of the spectrum of coronas lends itself to finding 
many large families of cospectral graph pairs. 



The published version of this paper appears in Linear Algebra and its Applications, Volume 435, 
no. 5, (2011), and contains occasional simplifications and additions beyond the current version. 
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Notation. The symbols 0„ and 1„ (resp., Omn and Imn) will stand for length-n 
column vectors (resp. m x n matrices) consisting entirely of O's and I's. For two 
matrices A and B, the matrix A^ B is the Kronecker (or tensor) product of A 
and B. For a graph G with adjacency matrix A, the characteristic polynomial of 
G is /g(A) := det(A/ — A). We use the standard notations P„, C„, Sn, and Kn 
respectively for the path, cycle, star, and complete graphs on n vertices. 

2. The Main Theorem 

Let G and H be finite simple graphs on m and n vertices, respectively, and 
let A and B denote their respective adjacency matrices. We begin by choosing a 
convenient labeling of the vertices oi Go H. Recall that G o H is comprised of the 
TO vertices of G, which we label arbitrarily using the symbols {1,2,..., to}, and 
TO copies Hi,H2, . . . , Hm of H. Choose an arbitrary ordering hi,h2, ■ ■ ■ ,hn of the 
vertices of H, and label the vertex in Hi corresponding to hk by the label i + mk. 
Below is a sample corona with the above labeling procedure: 






G 



H 



GoH 



Under this labelling the adjacency matrix of G o is given by 



AoB:= 



In® Im B®I„ 



The goal now is to compute the eigenvalues of this corona matrix in terms of the 
spectra of A and B. We introduce one new invariant for this purpose. 

Definition 1. Let _ff be a graph on n vertices, with adjacency matrix B. Note 
that, viewed as a matrix over the field of rational functions C(A), the characteristic 
matrix XI — B has determinant det(A7 — B) = ///(A) ^ 0, so is invertible. The 
coronal xh{X) € C(A) of H is defined to be the sum of the entries of the matrix 
(A/ — B)~^ . Note this can be calculated as 

Xff(A) = l^(A/„-i?)^il„. 

Our main theorem is that, beyond the spectra of G and H, only the coronal of 
H is needed to compute the spectrum of GoH. 

Theorem 2. Let G and H he graphs with m and n vertices. Let Xh{X) he the 
coronal of H. Then the characteristic polynomial of G o H is 

/Goff(A) = /ff(A)"/G(A~Xff(A)). 

In particular, the spectrum of Go H is completely determined by the characteristic 
polynomials fa and fn, and the coronal xh of H. 
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Proof. Let A and B denote the respective adjacency matrices of G and H. We 
compute the characteristic polynomial of the matrix Ao B. For this, we recall two 
elementary results from linear algebra on the multiplication of Kronecker products 
and determinants of block matrices: 

• In cases where each multiplication makes sense, we have 

M1M2 (g) M3M4 = (Ml (g) M3)(M2 (g) M4). 

• If M4 is invertible, then 

det (^JJ^ = det(M4) det(Mi - M^M^^M^). 

Combining these two facts, we have (as an equality of rational functions) 

/Goif(A) = det(A/™(„+i) - AoB) 

Xlmn - B ® 
(Xln -B)®I„ 



Xlm 


-A 


— In 




Xlm 


~ A 


-In 





= det 
= det 

= det((AJ„ - B)®lm) det 



(AJ„ -A)- {\l ® Im){{XIn -B)® Im)~^{ln ® Im) 

= det(AJ„ - Br dct(A/„ -A- (l^(A/„ - B)-^ln) Im) 
= det(A/„ - B)" dct((A - XH{X))Im - A) 

= fH{X)^fG{X-XH{X)). 

□ 

Remark 3. A natural question is whether or not the spectrum of G o is deter- 
mined by the spectra of G and H, i.e., whether knowledge of the coronal is strictly 
necessary. Indeed it is: Computing the coronals of the cospectral graphs ^5 and 

C4 U ifi , we have 

5A + 8 5A-2 

XS5(A) = 72 7 and XC4Uifi(A) 



A2_4 — A04u^iv'v A2_2A- 
Thiis cospectral graphs need not have the same coronal, and hence for a given graph 
G, the spectra of G o 6*5 and G o (G4 U Ki) are almost always distinct. Note that 
this stands in stark contrast to the situation for the Cartesian and tensor products 
of graphs. In both of these cases, the spectrum of the product is determined by the 
spectra of the components. 

The unexpected simplicity of the examples in the above remark are representative 
of a fairly common phenomena: since the coronal Xh{X) = ^f|xy ^^^^ computed 

as the quotient of the sum xh{X) of the cofactors oi XI — B by the characteristic 
polynomial /i/(A), it is a priori the quotient of a degree n— 1 polynomial by a degree 
n polynomial. In practice, however, as in the examples in the remark, these two 
polynomials typically have roots in common, providing for a reduced expression for 
the coronal. Let us suppose that ^(A) := gcd(xi/(A), ///(A)) has degree n — d (the 
gcd being taken in C[A]), so that Xh{X) in its reduced form is a quotient of a degree 
d — 1 polynomial by a degree d polynomial. Since the denominator of this reduced 
fraction is a factor of ///(A), and since fa is of degree m, each pole of xh{X) is 
simultaneously a multiplicity-m pole of /g(A — Xh{X)) and a multiplicity-m root of 
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fHi^)™" ■ Since these contributions cancel in the overall determination of the roots 
of /Goff(A) in the expression 

/GoH(A) = /H(Ar/G(A-XH(A)) 

from Theorem [2l we can now more explicitly describe the spectrum of the corona. 
Namely, let d be the degree of the denominator of X//(A) as a reduced fraction. 
Then the spectrum oi G o H consists of: 

• Some "old" eigenvalues, i.e., the roots of /ff (A) which are not poles of xh(A) 
(or equivalently, the roots of (7(A)), each with multiplicity |G|; and 

• Some "new" eigenvalues, i.e., the values A such that A — X-fr(A) is an eigen- 
value /i of G (with the multiplicity of A equal to the multiplicity of /i as an 
eigenvalue of G.) 

Since for a given /i, solving A — X_ff(A) = ^ by clearing denominators amounts to 
finding the roots of a degree d + 1 polynomial in A, the above two bullets combine 
to respectively provide all {n — d)m + m(d + 1) = m(n + 1) eigenvalues oi G o H . 

The following table, computed using SAGE ([H]), gives the number of graphs 
on n vertices whose coronal has a denominator of degree d (as a reduced fraction), 
as well as the average degree of this denominator, for 1 < n < 7. 



Table 1: Number of graphs on n vertices 
whose coronal has denominator of degree d 



d\n 


1 


2 


3 


4 


5 


6 


7 


1 


1 


2 


2 


4 


3 


8 


6 


2 







2 


5 


12 


28 


44 


3 









2 


13 


50 


138 


4 











6 


40 


304 


5 













22 


246 


6 












8 


214 


7 














92 


Total 


1 


2 


4 


11 


34 


156 


1044 


Average d 


1 


1 


1.5 


1.82 


2.65 


3.41 


4.68 


(Average d)/n 


1 


0.5 


0.5 


0.45 


0.53 


0.57 


0.66 



Since determining the characteristic polynomial oi G o H from the spectra of 
G and H requires only the extra knowledge of the coronal of H , it remains to 
develop techniques for computing these coronals. In Section 3, we will develop 
shortcuts for these computations, but we briefly conclude this section with some 
more computationally-oriented approaches. A first such option is to have a software 
package with linear algebra capabilities directly compute the inverse oi XI — B 
and sum its entries, as done in the computations for Table 1. This seems to be 
computationally feasible only for rather small graphs (e.g., n < 12). A second, 
more graph-theoretic, option relies on a combinatorial result of Schwenk j9j to 
compute each cofactor oi XI — B individually, before summing them to compute 
the coronal: 



SPECTRA OF CORONAE 



5 



Theorem 4 (Schwenk, 9 ). For vertices i and j of a graph H with adjacency 
matrix B, let Vij denote the set of paths from i to j . Then 

adj(A/ = ^ /h-p(A). 

Again, this approach becomes computationaUy infeasible fairly quickly without a 
method for pruning the number of cofactors to calculate. We explore this idea in 
the next section. Regardless, from Theorem HI we obtain: 

Corollary 5. The spectrum of the corona G o H is determined by the spectrum of 
G and the spectra of the proper subgraphs of H (or more economically, only from 
the spectra of those subgraphs obtained by deleting paths from H). 

3. Computing Coronals 

In this section, we will compute the coronals Xh{^) for several families of graphs, 
and hence for such H obtain the full spectrum of the corona G o H for any G. The 
principal technique exploits the regularity or near-regularity of a graph in order to 
greatly reduce the number of cofactor calculations (relative to those required by 
Theorem 2]) needed to compute the coronal. In particular, we use these ideas to 
compute the coronals of regular graphs, complete bipartite graphs, and paths. 

For graphs that are "nearly regular" in the sense that their degree sequences are 
almost constant, we can take advantage of linear-algebraic symmetries to compute 
the coronals. We begin with two concrete computations, those corresponding to 
regular and complete bipartite graphs, before extracting the underlying heuristic 
and applying it to the coronal of path graphs. The case of regular graphs, first 
addressed in [I], is particularly straight-forward from this viewpoint. 

Proposition 6 (Regular Graphs). Let H be r -regular on n vertices. Then 

A — r 

Thus for any graph G, the spectrum of G o H is precisely: 

• Every non-maximal eigenvalue of H , each with multiplicity \G\. 

• The two eigenvalues 



/X -I- r ± ■\/(r — nY + 4n 
2 

for each eigenvalue /i of G. 

Proof. Let B be the adjacency matrix of H. By regularity, we have = rln, 
and hence (A/ — i?)l„ = (A — r)l„. Cross-dividing and multiplying by 1^, 

A — 7' A — r 

The only pole of xh{X) is the maximal eigenvalue X = r oi H, and the "new" 
eigenvalues are obtained by solving A — — fi for each eigenvalue fj, of G. □ 

It is noteworthy that all r-regular graphs on n vertices have the same coronal, 
especially given that the cofactors of the corresponding matrices (S — A/)~^ appear 
to be markedly dissimilar. The simplicity of this scenario, and the easily checked 
observation that cospectral regular graphs must have the same regularity, lead to 
the following corollary. We will make use of this corollary in the final section. 
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Corollary 7. Cospectral regular graphs have the same coronal. 



As a second class of examples, we compute the coronals of complete bipartite and 
rt-partite graphs. 

Proposition 8 (Complete Bipartite Graphs). Let H = Kp^g be a complete bipartite 
graph on p + q — n vertices. Then 



XHiX) = 



nX + 2pq 



— pq 

For any graph G, the spectrum of G o H is given by: 

• The eigenvalue vuith multiplicity m{n — 2); and 

• For each eigenvalue pL of G, the roots of the polynomial 



/ix^ ~ {p + q +pq)x + — 2). 



Proof. Let B = 



0. 



0„ 



be the adjacency matrix of Kp q and let X = diag(((7 + 

A)/p, {p + \)Iq) be the diagonal matrix with the first p diagonal entries being {q + X) 
and the last q entries being (p + A). Then {XI — B)Xln = (A^ — pqjln, and so 



xh(A) = i^(a/-b)-4, 



l^Xln {p + q)X + 2pq 



A2 



pq 



A2 



pq 



Thus the coronal has poles at both of the non-zero eigenvalues ±^/pq of Kp^q, leaving 
only the eigenvalue with multiplicity p + q — 2. Finally, solving A — Xh{X) = /i 
gives the new eigenvalues in the spectrum as stated in the proposition. □ 

Remark 9. It might be tempting in light of Propositions [6] and [8] to hope that 
the degree sequence of a graph determines its coronal. This too, like the analogous 
conjecture stemming from cospectrality (Remark [3]), turns out to be false: The 
graphs P5 and K2 U K3 have the same degree sequence, but we find by direct 
computation that 



5X^ + 8A - 1 
A3 - 3A 



XK2UK3{X) 



5A-7 
A2 - 3A + 2 ■ 



Similar to the complete bipartite computation, we have the following somewhat 
technical generalization to complete /c-bipartite graphs. 



Proposition 10. Let H be the complete k-partite graph K„-^.n2,....nk- 
( 



Then 



XhW 



- 1 



A) 



V 



1=1 



A^-E-=2(j-i)QA'=-^ 



where Cj is the sum of the (^■) products of the form Ui-^i 
indices. 



■ Hi . with distinct 
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Proof. Let B be the adjacency matrix of H, let gj{\) = Yl i^i 

1=1 

the block diagonal matrix whose j-th block is gjluj ■ Then 



A) , and let X be 



k k 

(A/„ - B)X1„ = Y[{n, + A) - ^ n,gj 
_j=i j=i 

Solving as in the bipartite case, we find 



Xh{X) = lliXI ~ B)-'U = 



Ek 
i=i nj9j 



nti( 



•^)-E,=i'^j5j' 



which gives the result. 



□ 



The proof technique for the last two propositions generalizes to "nearly regular" 
graphs, by which we mean graphs H for which all but a small number of vertices 
have the same degree r. In this case, we can write 

(A/-B)l„ = (A-r)l„+v, 

where v = (vi) is a vector consisting mostly of O's. This gives 



(A/-B)-il„ 



1 



A 



[l„-(A/-S)-iv] 



and thus, using the adjugate formula for the determinant. 



xh(A) = i^(a/-b)-4, 



1 



A - r 



1 



fHiX) 



l<i,j<n 



where Cij denotes the (i, j)-cofactor of XI — B. Since Vi is zero for most values of 
i, we have an effective technique for computing coronals if we can compute a small 
number of cofactors (as opposed to, in particular, computing all of the cofactors 
and using Theorem 14]). For example, if we let /„ = /p„(A) be the characteristic 
polynomial of the path graph P„ on n vertices (by convention, set /o = 1), we can 
compute coronals of paths as follows: 

Proposition 11 (Path Graphs). Let H = Pn- Then 



XH 



(A - 2)/„ ■ 

Proof. In the notation of the discussion preceding the proposition, we take r = 2 
and V = [1 • • • 1]^. Further, we note that an easy induction argument using 
cofactor expansion gives Cja = Cj^n = fj-i- Thus we obtain 



XP„(A) = 



1 



1 



1 



A 



ri-l 



^" j=0 



A- 2 

from which the result follows. 

From this, we easily calculate the coronals for the first few path graphs: 



□ 
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Table 2: The coronals xp„(A) for 1 < n < 7 



n 


1 2 3 4 5 


6 


7 




1 2 3A+4 4A+2 5A^+8A-1 


6A^+4A-4 


7A^ + 12A-^-6A-8 


A A-1 A2-2 A2-A-1 A3-3A 


A3-A2-2A+1 


A'i-4A2+2 



Remark 12. This particular example can also be computed using Theorem For 
i and j distinct, there is a unique path from vertex i to vertex j, so the sum 
in the theorem reduces to a single term: 

adj(A/ - B)ij = fp^^[ij] = fi-ifn-j- 

Similarly, we find adj(A/ — B)ii — fi-ifn-i- Summing over all the cofactors gives 

if" " \ 

XP„(A) = - ^/.-i/„-. + 2 J2 , 

which reduces to the result computed in Proposition [11] after some arithmetic. 

4. COSPECTRALITY 

At the end of [1], the authors note that if Gi and G2 are cospectral graphs, then 
Gi o Ki and G2 o Ki are also cospectral, and that (by repeated coronation with Ki ) 
this leads to an infinite collection of cospectral pairs. Armed with the characteristic 
polynomial 

/Goif(A) = /H(Ar/G(A-XH(A)) 
of the corona (Theorem[2|), we can greatly generalize this observation on two fronts. 

Corollary 13. // Gi and G2 are cospectral, and H is any graph, then Gi o H and 
G2 o H are cospectral. Further, if Hi and H2 are cospectral and XHi = XH2 , o,nd 
G is any graph, then G o Hi and G o H2 are cospectral. 

We remark that examples of this second type do indeed exist. Define the switch- 
ing graph Sw(T) of a tree T with adjacency matrix At to be the graph with 
adjacency matrix 



^Sw(T) 



"1 


0' 


(E)At + 





1" 





1 


1 






and let Ti and T2 be non-isomorphic cospectral trees with cospectral complements 
(note that by [5], generalizing [5], "almost all" trees admit a cospectral pair with 
cospectral complement). Then the switching graphs Sw(ri) and Sw(r2) are non- 
isomorphic cospectral regular graphs (see fG*, Construction 3.7), and also have the 
same coronal by Corollary [Tj Corollary [13] now implies that for any graphs G and 
H, we have the cospectral pair G o Sw(T'i) and G o Sw(T2) and the cospectral pair 
Sw(Ti) o H and Sw(r2) o H. This gives, for example, infinitely many cospectral 
pairs of graphs with any given graph G as an induced subgraph. 
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